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Q | The reaction pd — > d<pp sp is studied within the Bethe-Salpeter formalism. 

^sO | Under special kinematical conditions (slow backward spectator proton p sp 

and fast forward deuteron) relevant for forthcoming experiments at COSY, 
I ■ the cross section and a set of polarization observables factorize in the contri- 
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for studying threshold- near processes at quasi- free neutrons. 
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3 I. INTRODUCTION 

a \ 

Data of elementary nucleon reactions involving a neutron in the entrance channel are 
^ scarce. One can use, for instance, a tagged neutron (n) beam [1] impinging on a proton (p) 

target or a proton beam impinging on a nuclear target. In the latter case one needs models 
to deconvolute the pn reaction. Here we are interested in threshold-near meson production. 
Over the last decade the data basis for the pp reactions has been impressively improved [2] . 
Such a detailed knowledge of elementary hadron reactions is necessary for analyzing heavy- 
ion collisions. Threshold-near meson production is also interesting for testing predictions of 
QCD in the non-perturbative region, provided e.g., by chiral perturbation theory. While it 
seems to be straightforward to translate cross sections for reactions of the type pp —>■ ppM 
(M is a meson) into the cross section nn nnM, the channel pn pnM is less simple 
accessible. Model calculations point to a non-trivial energy dependence of the ratio a(pp — > 
ppM)/a(np — * npM), e.g. for <fi meson production [3]. Experiments showed, furthermore, 
a surprisingly large ratio a(pd — > pdn)/a(pp — > ppn) [1,4]. Therefore, it is desirable to test 
directly predictions for the pn reactions. 

The spectator technique [5,6] offers one possibility to isolate to a large extent the quasi- 
free reaction at a neutron when a proton beam impinges on a deuteron target. Since the 
final state interactions in the pn system differs for singlet and triplet states and can be large 
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at threshold, it is of great advantage to consider the particular reaction pn — > dM, as the 
deuteron d is in one well defined state. Guided by this, in [7,8] the reactions pn — > dM with 
M — rj, n', u, 0, Oq have been considered. 

The aim of the present note is to analyze the complete reaction pd — > dMp sp and to 
test to what extent the spectator technique [5,6] really selects a quasi-free reaction at the 
neutron in the deuteron target. In other words, we are going to derive a factorization 
theorem showing kinematical conditions for which the total cross section a(pd — * dMp sp ) 
factorizes in a part depending on the target deuteron characteristics and a part a(jm — > dM). 
In doing so we select as a sufficiently transparent example the production. The reasons 
are obvious: (i) the elementary production amplitude is simple [3], (ii) the meson is 
interesting with respect to OZI rule [9] and hidden strangeness in the nucleon [3,8,10,11], 
and (iii) the meson is important with respect to both the inclusive K~ production [12] 
and the e + e~ decay channel to be studied at HADES in near future [13]. An extension to 
p and uj production goes along the same lines, but due to the role of s, u channel baryon 
resonances [10,14,15] the elementary amplitude of the subprocess pn — > pnM is by far more 
complicate and needs separate consideration, as other mesons too. 

Our paper is organized as follows. In section II we present the kinematics, cross section 
and amplitudes of the complete reaction pd — > d<pp sp . Numerical results are exhibited in 
section III, where polarization observables are discussed and compared with the ones of the 
quasi-free sub-reaction pn — ► d<p. We also derive selection rules. The factorization of the 
cross section is considered in section IV. The summary can be found in section V. 



II. KINEMATICS, CROSS SECTION, AMPLITUDES 

Let us consider the process of meson production in the exclusive reaction envisaged in 

[6] 

Pp+Pd = Pd + P4>+Pp, (!) 

at kinetic energies of the incoming proton corresponding to those achievable at the Cooler 
Synchotron COSY T p < 2.7 GeV. The deuteron is here supposed to be detected in the 
forward direction, i.e., in direction of the incoming proton with relativistic energy, T' d ~ 1 
GeV, and, contrarily, the final (spectator) proton is slowly moving in the backward direction 
(6 sp > 90°, p' p < 100 MeV). The main contribution to the process (1) comes then from the 
spectator mechanism, where the incoming proton interacts with the internal neutron from 
the deuteron and produces a meson and a deuteron in the final state. The second proton in 
the deuteron acts merely as a spectator. Schematically the spectator mechanism is depicted 
in Fig. 1. 

In the laboratory system, where the target deuteron is at rest, the relevant momenta are 
defined by p p = (E p , p ) and pd = (M d , 0) as initial momenta of the proton and deuteron, 
p' p = (E' p ,p' p ), p' d = (E' d ,p d ), and p^ = (E^p^) as momenta of final proton, deuteron and 
meson, respectively. The fivefold differential cross section reads 

OF' db(T - 1 
d d 3 p' d dQ' p ~ 8(2tt)VX 
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where the Kallen function A = 4 Mj \p p \ 2 and the kinematical coefficients B = Ma + E p — E' d 
and C = \p | cosO'p — \p' d \ cos9' d reflect the energy conservation. The invariant amplitude 
T depends, besides the spin projections s p , M.^ s' p , M.' dl M.^ on the quantization z-axis (ori- 
ented along p p ), also upon the initial energy and five independent kinematical variables, say 
p^ and Q! p . 

To compute the amplitude J 7 , let us define the scattering operator Om^,M' ( a 4®4 matrix 
in the spinor space) as the truncated amplitude for the subprocess pn — > d(p (hatched blob 
of the diagram in Fig. 1), i.e., as an operator which, when sandwiched between two free 
spinors u p and u n , results in the amplitude of the real process pn — > d<p with two on-mass 
shell nucleons. In the general case this operator is assumed to describe the process of meson 
production from two off-mass shell nucleons. With this definition, the amplitude T may be 
written as 

T = U a (p' p , S p ) (p' p - m) af) <S>M d (Pn,P p ) OM^M' d U *(PP' S p)' ( 3 ) 

where m is the nucleon mass and $^ d (j9 n , p' ) denotes the Bethe-Salpeter (BS) amplitude of 
the deuteron in the momentum space. In Eq. (3), the spinor indices are exhibited explicitly, 
with Latin and Greek characters referring to neutrons and protons, respectively (summation 
over indices occurring pairwise is supposed). An expression for the amplitude (3) can be 
found by specifying a model for the operator O. However, in order to keep the correspon- 
dence between the reaction (1) and the process of meson production in NN collisions as 
close as possible, we highlight explicitly its dependence upon the BS amplitude of the final 
deuteron by introducing another operator O as 

6%,m=-' J ^*%. d {l',r)d%(12;l'*,M+). (4) 

Then from Eq. (4) and from Fig. 1 it is easy to recognize the operator 0(12; 1'2', M.^) as the 
one describing the meson production in the sub-reaction pn — > d(p with an off-mass shell 
initial neutron. Such an operator has been introduced in Ref. [16] where its properties have 
been discussed in detail. We briefly recall here the main steps in computing 0(12; 1'2', M.^) 
and the amplitude T in Eq. (3). Usually, within the BS formalism in order to exhibit 
formulae in a familiar matrix form, one introduces new, "charge-conjugated" quantities 
such as * = -$7 C , Q = j c 6, -f c rf = - 7m7c , v(p) = [u{ph c } T , v(p) = [l c u{p)} T . Then the 
amplitude T can be written as 

T = v{p p , S p ) Q M ^M' d ^M d (Pn,Pp) (Pp + m) v{p' p , Sp), (5) 

where in Eq. (5) and in subsequent expressions one has usual matrix multiplications and 
traces. For further simplification we insert between Q and ^ in (5) the complete set of the 

4 

Dirac free spinors I = (l/2m) J2 £ r u r (p)u r (p), where e r — +1 for r = 1, 2 and e r = — 1 for 

r=l 

r = 3, 4, to obtain 

F= J-^2 T £?,t Md (v P n,t) u{p n ,s n )ty Md {p n ,p' p ) {p' p + m)v{p' p ,s' p ). (6) 
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The amplitude T^fj l Md (v pn ,t) = v(p p: s p ) Qm^M' u (Pn, s n ) coincides with the matrix ele- 
ment of the real process pn — > d<p, defined by the Mandelstam variables v pn = (p p +p n ) 2 , 
t = (p p —p' d ) 2 . Here is worth emphasizing that the kinematics of the process (1) differs from 
the one of the real process pn — > d<f>. In spite of commonly defined variable t, the numerical 
values of momenta of the final meson and deuteron do not coincide in the two processes. 
Rather, the amplitude T^ l ,t^ Md (i/ pn ,t) in the process (1) corresponds to an amplitude of 
meson production from an off-mass shell neutron. If the off-mass shell effects were negligi- 

'pn,t) in both processes would depend solely upon the momentum 



bly small, then T^ ) t^ Ai [y v 
transfer t and the effective energy vnn = (E p + E')' 

T^ti l Md can be written in the form 



p') 2 . As shown in Ref. [16], 
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r ,(12;l'2')^(2')^(l', 2')«r(l'), 
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where the partial spin amplitudes X ^.^,(12; 1'2') can be calculated within an effective 
one-boson exchange model by an explicit evaluation of the elementary diagrams depicted 
in Fig. 2 in [16]. We use the parameter set "B" in Ref. [11] to calculate the amplitudes 
A * n . rr t(12; 1'2'). Now the second part of the matrix element in Eq. (6), which includes the 
BS amplitude ^M d f° r the deuteron at rest, can be straightforwardly evaluated by using the 
explicit expressions for the 3 Si + and 3 Df + components of the deuteron BS amplitude [17]. 
The result is 



M*,M' dl 



y pn 



, (8) 



where M d is the deuteron mass, us,d denote S and D waves the BS wave function [17,18], 
n = P/|p| an d ^2m(n) the spherical harmonics of second order. In deriving Eq. (8) we 
used the identity \pl Ci„, i„„ Xi M Xi„„ = ^{ (T ^M d ) a 'y^ where C\_A.. are the Clebsch- 
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where C\ M t 

2"1 2^2 



Gordan coefficients, x\ v anc ^ ^ stand for the spin-| functions and the deuteron polarization 
vector, respectively. Equation (8) shows that an exact factorization of the cross section into 
two parts (one depending upon the target deuteron wave function, another one upon the 
kinematics and dynamics of the subprocess pn dip) generally does not occur. This is 
because the summation over the spins of the internal neutron s n which leads to interference 
terms in the square of the matrix element. However, below we show that under special 
kinematical conditions an approximate factorization is possible. 



III. NUMERICAL RESULTS 

From Eqs. (6-8), the differential cross section (2) and any other polarization observables 
of the process pd — > d<pp sp are calculable. The BS wave functions us,d basically coincide 
with the non-relativistic S and D waves of the deuteron for small values of the momentum 
of the spectator proton (see [17]). Therefore, the relativistic and non-relativistic expressions 
for the amplitude T have the same form, and one may perform also calculations by using 
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deuteron wave functions obtained within a non-relativistic approach. In our calculations we 
use the BS solution [18], obtained with a realistic one-boson exchange kernel. Additionally, 
Paris and Bonn wave functions are considered to investigate the sensitivity of observables 
on the deuteron model. 

All our calculations have been performed within kinematical conditions corresponding 
to the COSY Proposal #75 [6], i.e., the detected deuteron is supposed to be in the forward 
direction and the spectator proton in the backward hemisphere with 9' ~ 120° — 140°. The 
momenta of the deuteron are chosen in the range p' d ~ 1.9 - 3 GeV/c, which means that the 
spectator proton is detected with low velocity, i.e., p' p ~ - 0.04 GeV/c. These kinematic 
conditions imply that both the reaction (1) and the subprocess pn — > d<p occur essentially 
near the threshold. 

In Fig. 2, the differential cross section is presented as a function of the deuteron momen- 
tum for the specific choice of variables 9' d = 0°, and 9' = 130°. (We changed the kinematics 
with 9' d up to 5° and find a very smooth behavior of all observable.) The kinetic energy of 
the initial proton is T p = 2.69 GeV. The maximum of the cross section at p' d ~ 2.35 GeV/c 
corresponds to that kinematical condition where the momentum of the spectator proton is 
maximal, p' p ~ 0.035 GeV/c. Together with the large spectator angles this implies that 
at these kinematical conditions the excess energy in the pn — > d(p subprocess is minimal 
(As 1//2 ~ 4 MeV) and very near the threshold. Obviously, since at such kinematics the pro- 
duction amplitude is almost constant, the maximum in the differential cross section has a 
pure kinematical origin being entirely governed by the kinematical factor within the square 
brackets in Eq. (2). Fig. 2 highlights that different models for the deuteron wave function 
result in quantitatively different cross sections. This is because the integration over the 
internal momentum in the final deuteron covers the region of the minimum of the S wave 
where different models provide essentially different wave functions (see discussion in [16]). 

Let us now discuss polarization observables. Since the spin structure of the amplitude T 
in Eq. (8) is basically determined by the spin structure of the amplitude T Sp ,t^ d (iy pn , t) of the 
subprocess pn — > d<f), we proceed with an analysis of this amplitude. Bearing in mind that 
T Sp ,tn d {v P n, t) is a Lorentz invariant quantity one can, for the sake of simplicity, investigate 
its structure in the center of mass of the detected deuteron (near the threshold this system 
coincides with the center of mass system of two nucleons), provided the Lorentz boost is 

performed along the z-axis. In such a case the spin structure of T Sp ,t^ d {y pn , t) is the same as 
in the laboratory system (no additional Wigner rotations are needed). Near the threshold, 
the final <fi — d system is in a S orbital state corresponding to If = 0, = 1 + (Ad = 0, 2), 
Jj = + J = 0~, 1~,2~, where /, L and J n are the total isospin, total radial angular 
momentum, angular momentum and parity, respectively. Thus from symmetry constraints 
the allowed configurations in the initial state are 1, = 0, = 1, 3, 5 • • • and total spin Si = 0. 
The conservation law J; = Jf implies Lj = 1, so that Jf = \~ . Hence, near the threshold the 

main contribution into the amplitude T Sp ,tn d (v pn ,t) comes from initial spin configurations 
with s p = — s n and, since the two nucleons in the final state form a deuteron (with the total 
spin S = 1), the most probable transitions near the threshold are expected to be those with 
nucleon spin-flips in the corresponding vertices and with the deuteron projections M! d = ±1. 
Consequently, the projections of the meson are to be M.^ = =Fl. 

Coming back to the complete reaction one can expect that the prevailing contribution 
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into the full amplitude T comes from configurations with Aid — s' p = —s p , if (i) the effects 
of D waves in the target deuteron may be disregarded near the threshold, (ii) the spectator 
mechanism describes the reaction appropriately, and (hi) the factorization (see below) holds. 
Hence, near the threshold the polarization observables of the reaction (1) are predicted to 
behave in the same manner as the ones of the subprocess pn — > dtp. Therefore, the study of 
various polarization observables allows to check the factorization. 
Let us define the asymmetry of the reaction as 



A 



da(sp + s n = 1) + dajsp + s n = -1) - da(s p + s n = 0) 
da(s p + s n — 1) + da(s p + s n = -1) + da(s p + s n = 0) ' 



(9) 



where the spin of the intrinsic neutron is defined as s n = Aid — s' p . Then, from the above 
selection rules near the threshold, the asymmetry A is predicted to be negative and close 
to -1. In Fig. 3, results of a complete calculation of the asymmetry A are presented. The 
kinematical conditions are as in Fig. 2. It is seen that in the whole kinematical range of p' d 
the asymmetry is basically the same as in the subprocess pn — > dip. Moreover, as expected, 
at p' d near the maximum the asymmetry approaches -1. 

Similar results are obtained for other polarization observables. For instance, the tensor 
analyzing power defined as 



T, 



20 



1 da(M'd = 1) + da(M'd = -1) - 2da(M' d = 0) 
71 da(M' d = 1) + da(M d = "1) + da(M' d = 0) 



(10) 



is predicted to be l/y/2 from the selection rules. Indeed, the full evaluation of Eqs. (6 -8) 
yields a value close to that (see Fig. 4). 

Apart from the polarization observables which near the threshold are determined directly 
by the corresponding quantities for the subprocess pn — > d(p, one can define another set of 
observables with completely different behavior in the two reactions, i.e., in the subprocess 
pn — > d<p and in the reaction (1). For this purpose let us define observables at fixed initial 



deuteron and proton spin projections, e.g., at Aid = and s f 



i 

2' 



respectively. Near 



the threshold this means that the spin projection of the neutron will be s n = —s' p and by 
measuring the spectator polarization one can form different combinations which determine 
(allowed or forbidden) transitions for the subprocess pn — > d<p. Together with conservation 
rules for the complete reaction this admits a study of different observables which vanish in 
the sub-reaction pn — > d(f>. Correspondingly, let us define the quantity 



W z 



da(M4, = 1) - da{Ai <i> = -1) 



da(M^ = l) + da(M (j> = -1) 



s p = i;A4 d =0 



(11) 



If one measures this quantity for the positive spin projection of the spectator, then in the 
subprocess pn — * dtp one has super allowed transitions (s p + s n = 0, Ai^ = ±1, Ai' d = 
=1=1) and, since in (11) a summation over the final deuteron projections is assumed, W zz is 
predicted to approach zero at the threshold (see [16]). 

A different situation occurs if one considers W zz at s' p — — 1/2. In this case, at least 
for t = t m in{max) (he., the forward/backward direction for the subprocess pn — > dtp), the 
conservation of the spin projection in the complete reaction (1) strongly suppresses the final 
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states with = M' d = ±1 and M.<j> = — 1, M' d = 0, so that the only contributing final 
states remain these with M.^ = 1, Ai' d = 0,-1 and M.^ = — 1, M.' d = 1, which imply 
W zz ~ 1/3, in contrast to the vanishing W zz in pn — > d(p [16]. Figure 5 illustrates our full 
numerical evaluation for W zz at s' p = —1/2. 



IV. FACTORIZATION 



Let us discuss in more detail the behavior of the cross section Eq. (2) with the matrix 
element (8) near the threshold. As mentioned above, the summation over the spin pro- 
jection of the intrinsic neutron s n in Eq. (8) prevents an exact factorization of the cross 
section into the cross section of the subprocess pn — > d<p and a "target factor" . However, 
from Figs. 3 - 5 we conclude that, since the deviations of the exactly calculated observables 
from those determined by the selection rules of the subprocess pn — > d<p are small, a fac- 
torization should hold indeed with a good accuracy. For instance, the asymmetry in Fig. 3 

demonstrates the validity of the selection rule s p = —s n , i.e., the amplitude T Sp ,£n d (]s pn ,t) 
is proportional to 5 Sp ^ Sn . Similarly, Fig. 5 shows that the main contribution comes from 

transitions with \A4' d \ = 1, and from the above selection rules the amplitude T^^ Md (u pn ,t) 
is also proportional to Sm' ,-av Hence, near the threshold one may represent the amplitude 

as T Sp ,s n d {Vpni = <W d -A^ $s p -s n T SpyM ' d (u pni t), where T SpMd (v pn , t) is the threshold-near 
amplitude of the process pn — > d(p. Then in Eq. (8) the summation over s n disappears and, 
as conjectured above, the unpolarized cross section Eq. (2) factorizes. Substituting this am- 
plitude into Eq. (2), introducing the corresponding flux factors for the sub- reaction pn — >• d(f) 
and performing the summation over all spin projections, the unpolarized cross section can 
be cast in the form 

n F > d 5 ° . 2 M^n-4m 2 )| P ;| 2 , 7/ 2/| n / h \ da(pn^dcj>) 

2Ed d^dW p ~ \p\\b\ P ' p \-ce>\ Mp.I) + MlPpDJ Jt ' (12) 

where the differential cross section da{pn — > d(p)/dt depends only upon the transferred 
momentum t (determined by p p and p' d ) and the invariant mass v pn (determined by p p and 
p' p due to the balance E n = M d — E' p , p n = — Pp). From Eq. (12) it is seen that by studying 
experimentally the amplitude (8), near the threshold, one can obtain directly information 
about the subprocess pn — > d<p. Here it is worth reminding that the factorization of the 
cross section holds only if the contribution of the D wave is negligible, i.e., if in Eq. (8) one 
can put s n = —s p and arrives at 



Fny/&M d T£t'£<(u pn ,t) u s (\p' p \) C^i_ s - u D (\p' p \) L _ s C 1 ^ Y 2m (n) 



1 VI 

mfj, 



(13) 



Eqs. (8, 12) and (13) show that the smaller the values of the spectator momentum are 
the higher the accuracy of the factorization is. This has been confirmed numerically by 
comparing exact results based on Eqs. (2, 8) and approximate results based on Eqs. (12, 
13). In the whole kinematical region of interests the difference is less than 0.5%. 

One may generally conjecture that a structure analog to Eq. (8) gives rise to factorization 
under suitable kinematical conditions, hence providing a theoretical foundation of quasi-free 
reactions at the neutron in the deuteron. 
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V. SUMMARY 



In summary, within a covariant Bethe-Salpeter approach we study the relation of the 
tagged quasi-free reaction pn — > d<p and the complete process pd — > d(pp sp near the threshold. 
Guided by the very proximity of the numerical results for cross section and polarization 
observables based on the exact formulae and on the factorized ones, we find that, for special 
kinematical conditions, selection rules are operative immediately causing a factorization 
a(pd — > d(pp sp ) = F a(jm — > dip) where F depends on deuteron and kinematical quantities. 
This is equivalent to an omission of the contribution of the deuteron D wave at small values 
of the spectator momentum. 

As an extension of the present work, an investigation of the twice tagged reaction dd — > 
nnMp S pip S p2 can be envisaged to study the quasi-free meson production by the nn reaction 
proposed in [19]. 
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FIG. 1. Diagram of the process (1). 
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FIG. 2. Differential cross section of the reaction pd - 
momentum p' d at p p = 3.505 GeV/c, 6d' = 0°, p / = 130°. 
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FIG. 3. As Fig. 2 but for the asymmetry A. 
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FIG. 4. As Fig. 2 but for T 20 . 



11 




12 



